Conserved Quantities from the Equations of 

Motion 

(with applications to natural and gauge natural theories of gravitation) 

M. Ferraris* M. Francaviglia} M. Raiteri* 
Dipartimento di Matematica, Universita degli Studi di Torino, 
Via Carlo Alberto 10, 10123 Torino, Italy 



Abstract 

We present an alternative field theoretical approach to the definition of 
conserved quantities, based directly on the field equations content of a La- 
grangian theory (in the standard framework of the Calculus of Variations 
in jet bundles). The contraction of the Euler-Lagrange equations with 
Lie derivatives of the dynamical fields allows one to derive a variational 
Lagrangian for any given set of Lagrangian equations. A two steps algo- 
rithmical procedure can be thence applied to the variational Lagrangian 
in order to produce a general expression for the variation of all quantities 
which are (covariantly) conserved along the given dynamics. As a con- 
crete example we test this new formalism on Einstein's equations: well 
known and widely accepted formulae for the variation of the Hamiltonian 
and the variation of Energy for General Relativity are recovered. We also 
consider the Einstein-Cartan (Sciama-Kibble) theory in tetrad formal- 
ism and as a by-product we gain some new insight on the Kosmann lift 
in gauge natural theories, which arises when trying to restore naturality 
in a gauge natural variational Lagrangian. 

1 Introduction 

A number of physically reasonable geometric definitions of conserved quantities 
in field theories may be found in literature; in the recent past the issue to define 
conserved quantities for Lagrangian field theories has been in fact investigated by 
many authors on the basis of different formalisms. Just to mention a few of them 
we recall the Lagrangian method based on Noether's theorem [TT^ PTT71 RTTl f^lj] . 
the Hamiltonian approach and the symplectic methods 0] El El El El 
1441 1531 I63| . the Hamilton- Jacobi-based techniques |101 1131 1151 I16j . and the 
formulations which are directly based on field equations |S1 151)1 151) . These are 
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just some of the many references which are most relevant for our purposes, 
whereby other important literature is quoted. Since a satisfactory review is out 
of the scope of this paper we apologize for not being able to quote all authors 
and sources of information and we refer the reader to the quoted papers and 
references therein. 

In view of the fact that there exists in literature a widespread family of 
somewhat unrelated definitions it would be thence rather important to frame, 
as much as possible, all possible definitions under a unique and single method of 
construction. Very far from reaching such a fundamental goal, the scope of the 
present paper falls nevertheless into this line of thought. Our hope is to shed at 
least some new light on a fully covariant description of conserved quantities in 
gauge natural field theories. We choose the so-called gauge natural framework 
since it represents a rich geometrical structure which encompasses all La- 
grangian field theories relevant to fundamental physics (see also [251 1271 Pffll59| ): 
this modern formalism is in fact well suited to describe in a single mathematical 
context both natural theories (the most famous of which is General Relativity) 
as well as gauge theories together with their possible couplings with Bosonic 
and Fermionic matter. This setting becomes therefore a powerful tool to frame 
conserved quantities in a geometric setting which is general enough to include 
all relevant physical fields: gauge natural theories play in fact a privileged role 
to this purpose since they are defined, from the very beginning, on the basis 
of the covariance properties of the fields under the appropriate transformation 
groups. 

We point out that the gauge natural formalism is essentially a Lagrangian ap- 
proach. Nevertheless field equations are more fundamental than the Lagrangian 
itself in the description of dynamics. Field equations, in fact, rule out those 
field configurations which are physically admissible and dictate their dynamical 
evolution; accordingly we would prefer to focus our attention directly on the 
equations of motion. 

However, it has to be remarked that the solutions of field equations represent 
just a small part of the information that can be extracted out of a field theory. 
Other relevant information is encoded in the field symmetries; the knowledge 
on how fields are dragged along infinitesimal generators of symmetries is hence 
a further fundamental detail. 

Luckily enough, the two main ingredients which have to enter into a satisfac- 
tory definition of conserved quantities, namely the field equations content and 
the symmetry information, can be joined together into a unique structure. Field 
equations are indeed described via the Euler Lagrange morphism e(L) which 
turns out to be a differential form on (some suitable jet prolongation of) the 
configuration bundle of the theory. The Euler-Lagrange morphism can be con- 
tracted with the Lie derivatives £^y of the fields y with respect to infinitesimal 
generators £ of symmetries (which define vertical vector fields). The resulting 
object V := — < e(L)\£^y > turns out to be a horizontal form which can be 
interpreted as a new Lagrangian. From now on we shall refer to L' as the varia- 
tional Lagrangian (associated to L and £). The variational Lagrangian, for the 
information it encodes, is then in a good position to represent the fundamental 
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object out of which we derive conserved quantities: being L' a Lagrangian it can 
in fact be handled by means of the powerful tools of the Calculus of Variations 
in jet bundles. 

Our recipe to define the variation of conserved quantities will be thence 
developed in two steps which are both canonically and algorithmically well- 
defined at the jet bundle level. The first one is nothing but the first variational 
formula applied to the variational Lagrangian. The variation of L' , through a 
well-known integration by parts procedure, splits canonically into the Euler- 
Lagrange part (which vanishes because of symmetry properties) plus a pure 
divergence term Div T . The so-called Poincare-Cartan morphism T which (non 
uniquely but in a sense canonically) enters this divergence can be expanded as 
a linear combination of the coefficients of the vector field £ (generating the flow 
of symmetries) together with their covariant derivatives up to an appropriate 
(finite) order. The second step consists in implementing the so-called Spencer 
cohomology [HHIE2] through repeated integrations by parts with respect to these 
covariant derivatives. In this way we end up with a (to — 2) form U, called the 
potential, where to is the dimension of spacetime, which we choose to be a 
canonical representative, at the bundle level, suited to define the variations of 
conserved quantities. These latter quantities are indeed obtained by integrating 
the pull-back of U along solutions on the appropriate (to — 2) -dimensional 
domains of spacetime. We point out that the possibility of selecting a canonical 
representative for the potential IA is clearly an essential task, since different 
representatives may lead to different notions of conserved quantities. 

Notice that only the variation SQ of a conserved quantity, rather than the 
conserved quantity Q itself, is here defined. However, one can argue that SQ 
rather than Q is the truly fundamental object. Physically speaking, in fact, con- 
served quantities are not absolute, since only differences of physical observables 
are endowed with a direct physical meaning. Therefore one has to somehow 
fix a "reference point" (e.g. a background) to calculate them. However, in the 
absence of a linear structure on the space of fields there is no canonical choice 
for such a "reference point". Moreover, the formal on-shell integration of SQ 
depends on how we move along a curve in the space £ of solutions. Starting from 
a given solution ip , different deformations S correspond to different paths in £ 
passing through the same point ipo, so that SQ measures the variation of Q when 
we move from ifo to a nearby solution in a given direction. Roughly speaking 
the choice of such a direction corresponds to a control mode of the fields at the 
boundary which is based on physical grounds (e.g. micro-canonical or grand- 
canonical ensembles) and it is related to the asymptotic behaviour of specific 
solutions and to the specific physical observable one is aimed to determine (see, 
e.g., H EE! USED! 023 021 )■ We then believe that, inside a formalism which 
we want to be as general as possible, both physical and mathematical reasons 
suggest to us that only variations SQ should be reasonably considered. 

Finally, we point out another property that a viable definition of conserved 
quantities has to fulfill (and according to which we have developed the present 
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formalism). Conserved quantities are mathematically denned objects which are 
built out of fields with the purpose of extracting information about the physi- 
cal properties of all solutions; they represent indeed physical observables such 
as mass, energy, momentum, angular momentum or gauge charges. For this 
reason the mathematical definition of conserved quantities, whatever formalism 
we implement to obtain it, must eventually get rid of the mathematical struc- 
ture we started from and, in the meanwhile, it must not depend on the specific 
configuration variables we have selected to describe the solution itself. Let us 
better explain this "philosophical" viewpoint by considering as an example the 
theory of gravitation. It is well-known that there exist different Lagrangian 
formulations of gravity (e.g. purely metric |47| . metric affine |24j . purely affinc 
|23| . tetrad formulations Ashtekar variables 0, Chern-Simons formulation 
£Q and so on); these are described by profoundly different Lagrangians which 
in many cases are "equivalent" in the sense that they determine equivalent field 
equations even if they do not merely differ for the addition of divergence terms. 
From a geometric viewpoint these are really different theories since they are 
based on different configuration spaces and involve different dynamical fields. 
Nevertheless, all the aforementioned theories, under suitable regularity condi- 
tions, are equivalent on shell (i.e. they generate essentially identical or at least 
isomorphic spaces of solutions, possibly under appropriate gauge reductions). 
Basically, there exists a rule (which can be many-to-one) which maps a solu- 
tion of one theory into a solution of another. 1 For each different formulation 
of gravity we expect then a different setting for the definition of conserved 
quantities. Nevertheless it is physically desirable that all different definitions of 
conserved quantities turn eventually out to coincide on-shell; otherwise we may 
risk to obtain a paradoxical conceptual result. Indeed, if on-shell equivalent 
gravitational theories admit solutions that correspond to the same spacetime, 
all physically meaningful observables must be related to the spacetime itself and 
should not depend on which particular variable we have used to describe it. For 
example, the energy enclosed in a (bounded) region surrounding a black hole 
solution must depend solely on the solution and not on the specific configuration 
variable (e.g. metric, tetrad or connection) we initially choose to describe it. 

The formalism developed in this paper partially avoids such a possible para- 
dox. In fact, it is based on the variational Lagrangian which is obtained via 
contraction of the equations of motion with the Lie derivatives of fields. If any 
rule exists to unambiguously map the equations of motion and the Lie deriva- 
tives of one theory into the corresponding objects of another theory, then there 
exists, of course, an unambiguous correspondence between the variational La- 
grangians, too. Accordingly, also the conserved quantities of the first theory 
will be mapped into the conserved quantities of the second, so that, on-shell, 
they will eventually give rise in both approaches to the same numerical value 
for each conserved integral. 

However, we have to stress that there exist Lagrangian theories which ad- 

1 For example purely metric, metric affine and purely affine theories are all related by a 
generalized Legendre transformation and the quoted map among solutions is just given by the 
Legendre map 1541 
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mit different symmetry groups even if they are equivalent as far as their field 
equations content is concerned. For instance, purely metric, metric affine and 
purely affine theories are natural formulations of gravity and they all admit 
the group of spacetime diffeomorphisms as the "natural" symmetry group. On 
the other hand tetrad or Chern-Simons formulations of General Relativity are 
truly gauge natural formulations so that they admit a much larger group of 
symmetries, which, besides spacetime diffeomorphisms, encompasses also gauge 
transformations. 

In general, in gauge natural theories it happens that pure spacetime trans- 
formations cannot be globally isolated (in non trivial topologies) from gauge 
transformations. This fact reverberates then into the indeterminacy in defining 
the Lie derivative of gauge natural objects with respect to arbitrary spacetime 
vector fields (while such an indeterminacy does not occur for natural objects). 2 
Many lifts of the same vector field can be thence chosen: they originate different 
Lie derivatives and, accordingly, different variational Lagrangians and different 
conserved quantities. This fact raises a very interesting issue about symmetries 
(which we will face up in the last section) . Namely, albeit we have a correspon- 
dence between equations of motion we, a priori, lack a strict correspondence 
between symmetries. Nevertheless we claim that, a posteriori, also the latter 
correspondence can be gained in our approach and such a goal is simply ob- 
tained by imposing the equivalence of the variational Lagrangians. In this way 
a preferred lift of spacetime vector fields can be ruled out in gauge natural the- 
ories: it is the lift which restores the naturality in the gauge natural variational 
Lagrangian. 

The present paper is organized as follows. In section [21 in order to make 
the paper self-contained, we shortly review the geometric framework of natural 
and gauge natural theories. In section [3] we present the theoretical formulation 
of the "two steps procedure" which leads to the definition of the variation of 
conserved quantities. The result achieved here is tested for the purely metric 
formulation of General Relativity (in section^} and for the gauge natural tetrad 
formulation of General Relativity (in section 0). No discrepancies are observed 
between the two different formulations owing to the specific form of the varia- 
tional Lagrangians. The last section (section^ is finally devoted to the analysis 
of Einstein-Cartan theory in the gauge natural framework. This theory is a gen- 
eralization of Einstein's theory and, in vacuum, it becomes on-shell equivalent 
to it. We have chosen to investigate this theory since it clearly features the 
indeterminacy we mentioned above when defining conserved quantities in gauge 
natural theories with respect to spacetime vector fields. It will become however 
clear how this indeterminacy can be a posteriori eliminated by restoring natu- 
rality in the gauge natural variational Lagrangian. This procedure selects in the 

2 E.g., spinors, which are truly gauge natural objects cannot be classically Lie— dragged 
along generic spacetime vectors, but a classical notion of Lie derivative exists along Killing 
vectors, i.e. infinitesimal symmetries of the spacetime metric. A general definition of Lie 
derivative of spinors can be achieved only in the gauge natural formalism (see |27ll41ll5HI 'l. 
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specific example the so-called generalized Kosmann lift as the preferred lift to 
deal with. As a by-product we gain some new insight on the mathematical jus- 
tification of the Kosmann lift previously introduced in [27|^2[5!|| in the domain 
of gauge natural theories as a generalization of an ad hoc procedure introduced 
by Kosmann in |58| to Lie drag spinors. 

2 Natural and gauge natural theories 

Natural theories geometrically formalize the physical principle of general covari- 
ance. According to |351 156j we say that a field theory is natural when: 

A: the configuration bundle is natural; 

B: the Lagrangian describing the theory is natural. 

Item A means that spacetime diffcomorphisms can be functorially lifted to the 
configuration bundle Y (i.e. the space where the dynamical fields ip take their 
values). Roughly speaking, we know how fields transform under changes of 
coordinates in spacetime. As a consequence of naturality, for each spacetime 
vector field £ there exists a canonical lift £ on Y (which projects onto £). Using 
the natural lift £ it is then meaningful to consider Lie derivatives of fields with 
respect to spacetime vector fields by setting £^(p := £tip = Tip o(-(oj). 

Item B means that all spacetime transformations, once lifted on Y, are 
symmetries. Accordingly, each lift £ is an infinitesimal generator of symmetries. 

Obviously all theories based on diffeomorphism invariant Lagrangians which 
depend on tensor fields, tensor densities and/or linear connections are natural 
theories according to the given definition. 

In order to encompass into a unique geometric formalism theories admitting 
both diffeomorphism invariance as well as gauge symmetries one is led to in- 
troduce gauge natural theories. In gauge natural theories we assume that there 
exists a principal bundle (P, M, p;G), called the structure bundle, where all in- 
formation concerning symmetries is encoded; gauge natural theories admit the 
group Aut(P) of all automorphisms of the structure bundle as group of sym- 
metries (pure gauge symmetries correspond to vertical automorphisms). This 
amounts to say that a theory is gauge natural (see [23 E3 ESj for a deeper 
geometric insight) when: 

A: the configuration bundle is gauge natural; 

B: the Lagrangian describing the theory is gauge natural; 

C: a linear connection T and a principal connection A on the structure bundle 
can be built out of the dynamical fields. 

Item A is a pure geometric requirement: it means that the automorphisms of 
the structure bundle P functorially induce automorphisms on the configuration 
bundle Y. Accordingly, projectable vector fields Sp on P canonically induce 
vector fields S on 7. Item B is instead of dynamical nature since it implies 
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that all such induced vector fields H are infinitesimal generators of symmetries. 
Therefore we should be somehow able to associate conservation laws to each 
of them leading to quantities which are eventually physically interpretable as 
observables. In doing that we shall need the Lie derivatives of fields with respect 
to vector fields Hp, which are defined by setting 

£s P i P '■= £~±P = Tip o £ — H o ip (1) 

where £ denotes the projection of Hp onto spacetime M. Notice however that 
the group Diff(M) is not canonically embedded into Aut(P). We know how 
fields transform as a consequence of a transformation in P but we do not know 
how fields transform under change of coordinates in spacetime. In other words 
Lie derivatives with respect to spacetime vector field cannot be, at least a priori, 
defined for a generic £. 

Finally, item C is a technical requirement: the two dynamical connections F 
and A are the mathematical tools which are necessary to provide covariance at 
each step of the geometric formalism used to generate conserved quantities. 

Remark 2.1 We stress that despite in gauge natural theories there exists no 
natural way to define the action of diffeomorphisms on the dynamical fields there 
exist, however, many global (but not canonical) ways to lift spacetime vector 
fields up to the structure bundle. One of this, for example, is the horizontal 
lift defined through the dynamical connection A. Even if this lift apparently 
seems to be the most "natural" way to define the lift of vector fields, it has 
been nevertheless shown |3] that it does not lead to physically acceptable values 
for conserved quantities in physical applications so that eventually one has to 
resort to some other lift. In the applications we shall deal with, it will be shown 
that the generalized Kosmann lift is the most "natural" one; see section |SJ ■ 

Remark 2.2 We remark that, as far as the geometric formalism for conserved 
quantities is concerned, we are not so much interested in the Lagrangian but 
rather in the variational equations ensuing from it. Therefore we could have 
weakened item B by just requiring that the dynamical equations describing the 
physical system are (gauge) covariant. Nevertheless, starting from a set of co- 
variant equations of variational nature it is always possible, at least in principle, 
to build out a covariant family of Lagrangians depending on a (dynamical) back- 
ground; see jSH] (see also [j3 E] where covariant Lagrangians for Chern-Simons 
theories are exhibited). In the sequel we shall implicitly assume that any one of 
these Lagrangians has been already selected, even if we shall not be interested 
into its explicit form. M 

3 Conserved Quantities from the Equations of 
Motions 

Let us consider a gauge natural theory geometrically described through a gauge 
natural bundle (Y, M, tt) and dynamically defined by a fc -order Lagrangian L : 
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J k Y — ► A m (M). In terms of fibered coordinates (x^, y l ) on Y we locally have 
L = L(j k y)ds where L is the Lagrangian density, ds = dx 1 A . . . A dx m is the 
standard (local) volume form on M and (j k y) stands for (y l , y* . . . , J/L „ fc ), 
i.e. the set of partial derivatives of fields up to order k included. We shall denote 
by J k Y the k jet prolongation of Y, by V( J k Y) its vertical tangent bundle and 
by V*(J k Y) the vector bundle dual to V(J k Y). Notation and definitions in this 
section follow closely to which we refer the reader for a full treatment and 
further details. 

A variation is a vertical vector field X on the configuration bundle, which 
can be locally described as X = X l -g~ = Sy 1 -^—; physically speaking it de- 
scribes a one-parameter deformation of the dynamical fields. Accordingly, we 
can consider the variation SxL of the Lagrangian along the flow of (the pro- 
longation of) the vector field X. It is well-known (see [33] |5S1 |HS]) that each 
Lagrangian L induces a unique (global) morphism, called the Euler- Lagrange 
morphism 

e{L):J 2k Y — ► V*(Y) ® A m (M) (2) 

together with a family of (global) morphisms (which depend on the Lagrangian 
and possibly on a connection L on M) called Poincare-Cartan morphisms 

T{L, L) : J 2k - X Y — >V*{J h - x Y)® k m ~ x {M) (3) 

The Euler-Lagrange morphism and the Poincare-Cartan morphisms are in fact 
defined so that the so-called first-variation formula holds for any deformation 
X on Y: 

S X L = < e(L) | X > + Div < JF(L, L) | j^X > (4) 

where < | > denotes the canonical pairing between differential forms and vector 
fields (in our case between elements of V*( J h Y) and elements of V( J h Y)). The 
formal divergence operator on forms is defined by 

Div(/)o/+V = d(/°jV), f:J h Y^K{M) (5) 

d{ ) being the exterior differential operator on forms and A(M) = ©fcA (M) 
denoting the bundle of forms over M. The Poincare-Cartan morphisms are 
uniquely defined only for k = or k = 1; for k = 2 they are not unique but 
still there is a canonical choice which is independent on any connection T. For 
k > 2 the Poincare-Cartan morphisms strongly depend on the choice of a linear 
connection; see [SJ. Nevertheless, in gauge natural theories the existence of 
a dynamical linear connection T is axiomatically required (see item C of the 
definition); in this case we shall implicitly assume that the Poincare-Cartan 
morphism entering the first variation formula (for k > 2) is the one induced 
by this preferred connection L. For this reason, from now on we shall omit to 
indicate in the notation the dependence on T of T{L,T) and we shall simply 
write F{V) . 

The Euler-Lagrange morphism (J2J , which can can be locally written as 

e{L)=e i {j 2k y)dy i ®ds (6) 
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encodes the information relative to the equations of motion for the dynamical 
fields. A critical section (or a solution) is a section ip : M — >Y, locally de- 
scribed as (p : x^ i — ► (x^, y l = ip l (x)), the prolongation of which belongs to the 
kernel of the Euler-Lagrange morphism, i.e.: 

e(L)oj 2k ip = =^ e k (ip i ,d li (p i ,...,d lil ... lialt (p i )=0 (7) 

Now, let us denote by Ep a projectable vector field on the relevant principal 
bundle (P, A/, p) of the theory. It locally reads as Sp = £ M d p + £ A pa (where 
PA ia a local basis for right invariant vertical vector fields on P). It canoni- 
cally induces a vector field E on the configuration bundle which is, by the very 
definition of gauge natural theory (see the previous section), an infinitesimal 
generator of symmetries. Through the vector field E we can define the (formal) 
Lie derivative £~,y of the fields. Since £^y : J X Y — ► V"(Y) takes value into the 
vertical bundle of Y it is meaningful to consider the contraction 

L'(L, E) = -< e(L) \£~y>: J 2k Y — ► A m (M) (8) 

which defines a horizontal form on the bundle J 2k Y which we shall call, from 
now on, the variational Lagrangian. 

Since we are assuming the configuration bundle to be a gauge natural bundle 
the Lie derivative £&y entering into the definition of the variational Lagrangian 
can be written as a linear combination of symmetrized covariant derivatives with 
respect to the dynamical connections (T,A); see [2H1EH1- Thereby expression 
ijBJ) can be locally written as 

L'(L,E) = {W^ + W^V Pi e + --- + W^V (pi ... PT) e (9) 
+W A £ A + W£V Pl t A + ... + w£- p ° V (pi ... Ps) ^} ds 

where W^, W? 1 , . . ., W£ 1 " p ', W A , W£, . . ., W£'" Pr are tensor densities with 
respect to automorphisms of the structure bundle (the pair (r, s) is called the 
order of the gauge natural bundle; see |56|): here and in the sequel round 
brackets around indices denote symmetrization. 

Whenever we have such a linear combination we can perform covariant inte- 
gration by parts to obtain for L'(L, S) an equivalent canonical expansion under 
the form: 

L'(L,E) = B{L,~) + Div S(L,E) (10) 

where the quantity £ (L, 5) is usually called the reduced current and vanishes on- 
shell, while the quantity B(L, S) is linear in S and it turns out to be identically 
vanishing along any section (see I38| ) . The identities 

B(L,E) = (11) 

arc called generalized Bianchi identities (indeed they reduce to the usual Bianchi 
identities in General Relativity and in gauge theories). Hence equation gives 
rise to a conservation law: 

Div£(L,S) = L'(L,E) ~ (12) 
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where ~ denotes equality on-shell. Nevertheless it is clear that Q121 is a trivial 
conservation law owing to the fact that the reduced current £ , which is built 
out of the coefficients (0 of the Euler-Lagrange morphism together with their 
covariant derivatives, is, as we said, vanishing on-shell. Hence equation l|12fl 
is of no utility by itself to describe physical properties of any given solution. 
Nevertheless, the variational Lagrangian becomes the starting point to define 
algorithmically the quantities we are interested in. The formalism basically 
develops in two steps. 

First step: integration by parts with respect to X. 

Let us consider a vertical vector field X on the configuration bundle, locally 
given as X — X l -J~ — 5y l -^~. By treating (jHJ as a new Lagrangian we can 
consider the variation SxL' and, accordingly, we can make again use of the 
first variational formula 10} . Notice that L' in general depends on the fields y % 
together with their derivatives up to some order h < 2k, but it also depends 
on the variables and £ A together with their derivatives up to the orders r 
and s, respectively; see @. 3 In the variation SxL' the components £ M and £ A 
are independent on fields and thereby are kept fixed, i.e. Sx^ = $x£, A = 0. 
In the sequel we shall also be interested in the case in which the components 
£, A , via some suitable "lift", are built out of the dynamical fields y % and the 
coefficients £ M together with their derivatives up to some fixed order, i.e. £ A = 
£, A (j a Di3 b £,^)- When such a dependence is inserted back into (JHJ) we end up 
with a variational Lagrangian which depends only on the dynamical fields and 
the components Thereby the whole formalism we are going to develop will 
simply work by setting £ A = 0. 

The first variation formula applied to the Lagrangian L' reads now as follows: 

5 X L' = <e(L')\X>+DivF(L',X) 

= DivT(L\X) (13) 

where T(IJ ', X) stands for < T{L')\j p X > for some suitable p. The latter equal- 
ity in l|13(l is due to the fact that, being the Lagrangian L' a total divergence (see 
@H0 ^s associated Euler-Lagrange morphism < e(L') \ X >—< e(Div£ ) | X > 
vanishes identically. Moreover the variation of (|12l) turns out to be a pure 
divergence, i.e.: 

S X L' = Div6 x £ (L,E) (14) 

Hence, comparing (|13|) with i|14|) we obtain a strong conservation law (i.e. a 
conservation law which holds true off-shell): 

Dw{f(L',X)-5 x £(L,E)\ =0 (15) 

3 The inequality h < 2k is due to the fact that a Lagrangian of order k can give rise to field 
equations of order lower than 2k; see, e.g. General Relativity or Lovelock metric theories 1181 
where k = 2 while field equations are second order only. 
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In particular, if the variation X is a solution of the linearized field equations, 
i.e. it is tangent to the space of solutions, then Sx£(L, 3) = and we obtain a 
weak conservation law: 

DivF(L',X)~0 (16) 



Second step: integration by parts with respect to 3. 

We shall now show that the (m — 1) form !F(L',X) — 5x£(L,E) is not only 
closed but it is also exact (off-shell) and we shall explicitly exhibit a global 
representative for its potential. This is a result which does not depend on the 
topology of spacetime and it does not depend on any given solution. Indeed 
all calculations will be performed on some jet prolongation of the configuration 
bundle Y. Only at the end of the calculations we shall pull-back the results 
along sections of the appropriate jet bundle obtaining in this way differential 
forms on the base manifold. Moreover the global character of the results we shall 
obtain will follow directly from the fact that the algorithm developed fulfills the 
covariance property at each step of its construction. 

First of all let us consider the morphism T(L' , X). Since L' depends on j h y, 
and j s £, A , with h < 2k, and the expression in the right hand side of 11311 
is obtained, in general, through h integration by parts, the quantity T{V , X) 
depends linearly on the independent components £ A together with their 
derivatives up to the orders, respectively, h + r — 1 and h + s — 1. Therefore it 
can be written as 

F(L',X) = {^+^ V «^ + ---+- F r"' / '* +P " lV (Pl...PVHr-0^ 

+F p A z A +rrv P1 e + ■ ■ ■+K pi - Ph+s - i ^( P ,.. Ph+s - 1 ^ A } d Sp 

where ds p = d p \ds. The coefficients in (|17H are tensor densities which depend 
on the dynamical fields y l together with their variations X 1 = 5y l up to some 
(finite) order. Notice that covariant derivatives can be defined with respect 
to the dynamical connections (see item C of the definition). Notice also that 
the coefficients j^P pi - p p an d j:PPi—Pi are Q f course symmetric in the indices pi 
but not with respect to the whole set of upper indices. However, whenever we 
have a linear combination of the kind (|17|l we can perform covariant integration 
by parts to obtain for the same quantity an equivalent linear expansion the 
coefficients of which are all symmetric with respect to upper indices, while the 
integrated terms are all pushed into a formal divergence. In this way expression 
<|17|) can be recasted as follows: 4 

F{L',X) = F(L',X)+Divll(L',X) (20) 

4 For the sake of completeness of this paper we give the explicit formula for the decompo- 
sition in the case h + r — 1 = 2, s = and r a M „ symmetric. We have in this case: 

H = ^ + fS - V M^' 31 + fl^ ] ) 

= 4 ap) + /< Q/3) (18) 

fj, & 
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= [f p e + fjf^v P1 e + ■■■+ ^■■■ pi+r - i) v (P ,., t+r _ 1 )^ 

+P A S A + f%"\^ A + ■■■+ ^A PP1 - Ph+s - l} ^( Pl ... Ph+ ^ A } ds P 
+d a U lpa] ds p 

Here and in the sequel square brackets around indices denote skew-symmetrization. 
(For the general setting of the above decomposition for any pair (r, s) we refer 
the reader to [HHD- We only stress that this kind of integration by parts is 
a well-defined global operation in the bundle framework. Globality is indeed 
assured by the direct use of covariant derivatives. Uniqueness of T instead 
descends from its symmetry properties and from the fact that Div is a nilpo- 
tent operator: Div o Div = 0. Therefore, it can easily checked that the term 
jF(L',X) in the right hand side of (|2C)|> coincides with the variation 5x£ 
of the reduced current and it is hence vanishing on shell. Indeed, inserting l|2U[l 
into H15[) and taking Div o Div = into account we have: 

Div{T(L',X)-6 x £(L,Z.)} =0 (21) 

Since the components £ M and £ A are arbitrary and since both terms T{L',X) 
and dx£(L,E.) are, by construction, symmetric in the upper indices, it is easy 
to verify that implies f(L',X) = S X £{L,E). 

Also the potential U features an expression of the kind: 

+ U [ rk A + ■ ■ • +4 Pffl(Pl -' Ph+ ^ 2) V( P ,.. Ph+a „ 2) e A (22) 

which is clearly not unique but only defined modulo closed forms. Nevertheless 
(see I38) 1 ) there exists a unique representative the coefficients of which have the 
maximal symmetry property: U^ papi ^" Pp = 0. From now on we shall select this 
representative which is, from a mathematical viewpoint, the canonical one. The 
physical viability of this choice will be tested in applications. 

Let us summarize. Starting from the variational Lagrangian (JHJ we have ob- 
tained a current: 

J r (L',X)=S x £(L,E)+T>ivU(L',X) (23) 

which is conserved on-shell (see (|16f) 'l and it is also exact on-shell if X is a 
solution of the linearized field equations, since in this case the term Sx£(L, 3) 

with = 1/34 CT/3]Q Rl a0 , ff = -4/3V ff 4 aCT]fl and 

u *p = _ 2/3v CT ^ir /3]CT } e + 4/3^ [ r /3]CT v tT e i (19) 

The generalization to s ^ is easily obtained replacing lower Greek indices with Latin indices, 
replacing the Riemann tensor with the field strength F and enlarging the covariant derivatives 
to act on both internal and spacetime indices. 
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in (|23[) is vanishing (notice, instead, that the current !F(L',X) — Sx£{L,E) is 
exact also off-shell). 

Let us now consider a section ip : M — ► Y . Its prolongation jcp to a suitable 
(finite) order can be used to pull-back expression 1|23[) onto the base manifold. 
We shall denote by !F(X,(p), £(E,tp) and U(X,ip) the pull-backs on M of the 
relevant quantities in i|23|) (from now on the indication L 1 between the round 
brackets will be omitted for the sake of simplicity). Given a hypersurface D 
in spacetime we define the variation 5xQd(z;<p) of the conserved quantity 
Qd(E, ip), relative to the set (S, (p., D) as follows: 

8 x Q D (E,<p) := f F(X,cp) = f S x £(E,p) + f U(X,<p) (24) 

JD JD JdD 

~ I U(X,p) (25) 

J 3D 

where the last equality holds true whenever <p> is a solution of field equations 
and X is tangent to the space of the solutions, i.e. it describes a one-parameter 
deformation of ip along nearby solutions. In particular if the region D is a 
(portion of a) Cauchy surface and the vector field £ is transverse to D we identify 
^Qd(2,^) in 1|24[1 with the variation of the Hamiltonian Sx'Hd^i ¥>)■ We 
stress that, for the time being, this is just a definition. Nevertheless the name 
Hamiltonian will be justified by the applications we shall explicitly consider. 
Indeed it will turn out that SxTi-D^T f) physically describes the evolution of 
the fields along the flow of the vector £. 

Remark 3.1 We remark that, even though we have assumed a gauge natural 
Lagrangian as the starting point for our framework, the whole theory is in fact 
developed from the Euler-Lagrange morphism JSJ . The Lagrangian is necessary, 
a priori, only in defining the dynamics via its equations of motion. This implies 
that it is of no interest whatsoever Lagrangian representative we choose in the 
cohomology class [L] the elements of which differ from each other only by the 
addition of divergences. Indeed all the representatives inside the class [L] give 
rise to the same equations. Even more, we can relax the hypothesis that the 
Lagrangian be gauge natural since only the equations of motion are required to 
be gauge invariant as well as diffeomorphism invariant. Hence, our formalism 
allows to encompass also more general classes of field theories, e.g. Chern- 
Simons theories. M 

Remark 3.2 Notice that, through a direct application of the Noether uheo- 
rem, starting from the Lagrangian L and from the infinitesimal generator of 
symmetries S, one is able to algorithmically construct a conserved Noether cur- 
rent £ (L, S) and to exhibit a superpotential V(L, S) such that the following two 
identities hold: 

Div£(L,E) =L'(L,E) ~0 (26) 
£(L,E)=£(L,Z) + DivV(L,E) (27) 
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where £(L,E) is the same reduced current appearing in lt25|) (see [2H f° r 
details). By taking the variation of 1)261) and taking 1)131) into account we then 
get: 

Div S X £{L, 3) = DivF(L', X) (28) 

so that 5 X £(L, S) and T(JJ , X) differ from each other for a form which (at least 
locally) is exact: 

T{L',X) = 5 X £ (L,E) + DwC{L,E,X) (29) 

IP 6 X 8(L, 2) + Div {S X V(L, 2) + C(L, 3, X)} 

Notice now that the divergence DivC cannot be vanishing in general. Indeed, if 
DivC vanishes, a comparison between (|29|) and i)23|) leads to the identification 
5 X V(L,E) = U(L',X), which cannot be true. We can infer this without any 
calculation. Indeed, it is well known that if we consider a Lagrangian L = 
L + Div a which differs from a given Lagrangian L only for the addition of a 
divergence, then the reduced current remains the same, i.e. £(Z,S) = £(L,E) 
while the Noether superpotential transforms as follows (see |40l 148) 1 : 

V(L,E)=V(L,E)+i £ a (30) 

On the other hand, being equation i)23)) built out directly from the equations of 
motion @ , is not affected by the addition of divergence terms to the Lagrangian. 
For this reason we see that the divergence term DivC(L, 5, X) cannot be vanish- 
ing in general and it must be sensitive to the representative we choose inside the 
cohomology class [L] a given Lagrangian belongs to, in order to counterbalance 
the transformation rule 1)30)1 into 1)29)1 . 

Roughly speaking, we could say that the formalism developed up to now, 
namely the construction of T)ivU(L' , X) in 1)23)) . is nothing but a recipe to algo- 
rithmically define the additional term DivC(L,5, X) which has to be added to 
the Noether superpotential in order to make the variation of conserved quanti- 
ties insensitive to the choice of a representative inside [L] . 

When dealing with applications we shall see that the fibered morphisms 
U(L',X) and S X V(L, 3) differ indeed for a term which is nothing but the co- 
variant Regge-Teitelboim boundary correction term; see |28l 1361 1401 1501 . 



4 Purely metric formulation of gravity 

Let us now consider the Hilbert Lagrangian density 

Lh = ^g^R^ifg) (31) 

(k = 87rG/c 4 ) or any other Lagrangian differing from it by the addition of 
boundary terms (e.g. the first order Einstein Lagrangian |22) . the first order 
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covariant Lagrangian |29U31ll3l)ll51| or the York's Trace-K Lagrangian |16U66p . 
Indeed, what is really relevant is only the expression of Einstein's field equations: 

_yJ G ^ = _Yj i R ^_l g ^ R \ = (32) 
from which we define the variational Lagrangian Q: 

L'(0 = ^ £ i9fW ds=^- ds (33) 

where £ = denotes a vector field in spacetime. 



First step. Let us consider a vertical vector field X = X^-g^— = 5gp V -^— on 
the configuration bundle Y — Lor(M) of Lorentzian metrics and let us perform 
the variation 5xL'(^). By taking into account the relations: 

^=V p (5up, u p MV ='y$ u ,-6lX )a (34) 

h% = \g P °(-Va6g^ + V^g va + V v Sg ail ) (35) 

(where denotes the Levi-Civita connection of the metric g) we obtain: 

S x L'(0 = & r F'V 1> g,j 1 X,j 2 Z)d8 (36) 

where: 

F> = ^{j 2 g,X)e +^ p X {j 1 g,j 1 X)Vxe +^\g,X)V {x<r) e (37) 
with: 



777 
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(38) 



Sg a p (/ [ V h ^+5 p V [7 ^ ] 
+ l9 aP (9 plX 62 ] +9 llX S$)) (39) 



2~ 

Second step. Integrating the expression Ij37(l by parts according to the formu- 
lae (fT%|l and ijT§|) we obtain 

F> = S x S 1 (L,0 + d p U^ (40) 

where: 

f(L,i) = ^-GlC (41) 
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Notice that the reduced current £(L,£) is clearly vanishing on-shell. Moreover 
the first term in the right hand side of l|42|l is nothing but the variation of the 
Noether superpotential (i.e., in our case, the variation of the Komar superpoten- 
tial; [57|). On the contrary the second term in the right hand side of (|42|l is the 
covariant Regge-Teitelboim correction term |1UI 1361 HP) , namely it is the term 
which has to be added to the Hamiltonian to define a well-posed variational 
principle. To understand this latter statement, let us analyse (|40() in details. 
Let us choose a (local) section g : x i— > g^{x) of the configuration bundle and 
let us make use of it to pull-back formula (|4U[) onto spacetime M. Let D be 
a (portion of a) Cauchy surface in M and let the vector field £ be transverse 
to D. According with our previous definition we identify the integral (of the 
pull-back) of l|40|) with the variation of the Hamiltonian: 



SxH D (£,g) = IJ^GZe\ds^ (43) 



i)D 



2k q 



Since the vector field £ is transverse to D we can drag (at least locally) the 
surface D along the flow of £, defining in this way a m-dimensional region of 
spacetime which is foliated into hypersurfaces diffeomorphic to D. By making 
use of the (3+1) ADM formalism we can rewrite i|43|) in terms of quantities 
which are adapted to the foliation (the details of such a calculation can be 
found in PIBIIS1EDJ). We obtain: 

SxH D {^g) = f {5NH + 6N a H a + [h aP \5P af) - [P al3 ]8h aP } d 3 x (44) 

JD 

where N and N a are, respectively, the lapse and the shift of the vector field 
£ = dt, TL and Ti. a are the usual Hamiltonian constraints of General Relativity, 
P a/3 is the momentum conjugated to the metric h a p on the hypersurface D 
and [h a p] and [P a/3 ] denote the right hand side of the Hamilton equations: 
£(h af} = [h a p] and £ i P a(3 = [P^] (see (TUJ EDJ)- We point out that the 
boundary terms arising in the variation of the first term in the right hand side 
of H43|) exactly cancel out the second and the third boundary terms. Hence 
we end up with the "pure" bulk term l|44|) which is the correct expression one 
would expect for the variation of the Hamiltonian for General Relativity; see 
|1(JI II II IT51 141)1 163| . This fact justifies, a posteriori, the definition of variation of 
the Hamiltonian we have previously attributed to formula l|43|) . 

We just outline that the variation of the energy 6xE]j(^, g) enclosed in the 
quasilocal region bounded by dD and relative to the pair (£,g) can be defined 
as the on-shell value of l|43|l . namely: 

6 x E D {i,g) = jf j* +^9^Su^e ] ]d Spl (45) 

where, we recall, for variation we mean the infinitesimal change of energy when 
we move along nearby solutions (different vectors X corresponding to different 
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paths in the space of solutions all passing through the same "point" g{x)). The 
(3+1) splitting of formula l|45[) can be found in 0U| an d for this reason we do 
not exhibit it here. In |4l 131 II 91 l3*Tl 141)1151]) the relations interplaying between the 
integrability conditions of l|45|l and boundary conditions on the metric and its 
derivatives are analysed in detail. We just point out that expression l|15|) with 
Dirichlet boundary conditions gives rise to the Brown- York quasilocal energy 
1^02]. Moreover the implications of (|45|l in relation with the first law of black 
hole thermodynamics have been discussed in HHI 1201 HO E01 E01 EHl It was 
there also emphasized how this formula sheds some light on the fundamental 
contribution of geometry in characterizing the entropy of black hole solutions 
and even of more general solutions of Einstein's equations. 

For all these reasons we consider expression lj4"5|) as an acceptable and viable 
definition for the (variation of) energy. 

5 Gravity in tetrad formalism 

Let us now consider the formulation of gravity in the tetrad formalism. Even if 
the tetrad formulation is in many cases equivalent to the purely metric formal- 
ism, it becomes necessary when dealing with spinor matter, so that it deserves 
an investigation in its own. 

It is well known that a good mathemathical arena to globally describe tetrad 
gravity is the gauge natural bundle framework (see ^fl ) . In this context 
the structure bundle of the theory is a principal bundle (P, M,p; SO(l, 3)) over 
spacetime M. 5 According to [23 a tetrad field is defined to be a section of a 
GL(4,R} bundle S which is the bundle associated to the bundle P x L(M), 
where L{M) denotes the frame bundle, via the left action 

p : (SO(l, 3) x GL(4, R)) x GL(4, R) — ► GL(4, R) 

p : (A, J;X) t-> A- X ■ ( J -1 ) (46) 

Fibered coordinates on T, are denoted by (x M , e^), where i = 0, . . .3. A tetrad 
field is then a map x i— > (x^e* — 9 l ^(x)). Given a tetrad field we can define 
a metric over spacetime through the rule g^ — TjijO'lOl, rj — diag(—l, 1, 1, 1), 
so that, a posteriori, the structure bundle (P, M,p; SO(l, 3)) can be identified 
with the bundle of orthonormal frames SO(M, g) which is the subbundle of the 
frame bundle L(M) formed by all ^-orthonormal frames. We stress that, in this 
framework, the index i of the tetrad is not merely a label denoting a set of four 
1-forms in spacetime. In this case tetrad fields would be merely a local basis of 
the tangent space (globality being recovered only for parallelizable manifolds) 
and the theory would just be a natural theory. Global frames, i.e. tetrads, are 
instead truly gauge natural objects: they are sensitive to the transformations of 
the structure bundle P. Each authomorphism of the principal bundle functo- 
rially induces, via the left action p, a transformation law acting on the tetrad. 

5 Since we shall not consider here the explicit coupling between gravity and spinor matter 
the structure group is assumed to be the pseudo orthogonal group instead of its twofold 
covering group Spin(l, 3) 
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This transformation law allows to canonically define the Lie derivative of the 
tetrad fields with respect to any vector field on the structure bundle. Indeed a 
projectable vector field Sp in the principal bundle, locally described as 

~p = e (x)a ll + e j {x)Pn , e j = t m m 

(having denoted by p = gd/dg, g 6 50(1,3) a basis for right invariant vector 
fields on P in a trivialization (x, g) of P) canonically induces the projectable 
vector field 5 on S given by 5 = £ M <9 M + £*■ e^d / 'de*. Through the vector field H 
we can define the formal Lie derivative ,£~e : J 1 !] — >V(E). It is intrinsically 
defined by the rule £~.e o j l := £^8 = TO o £ — 5 o 6, where £ denotes the 
projection of S onto M and 9 is a tetrad field. In coordinates: 

£E^^Cd^ + d^el~e^ (48) 

In the tetrad affine formulation of gravity the configuration bundle Y is assumed 
to be the bundle: Y = E — > M. The Lagrangian of the theory turns out 
to be the fibered morphism L : J 2 E — >A 4 (Af) locally described through the 
Lagrangian density: 

^ = ^eiR k a l (j 2 e)e^e ljkl (49) 

where R kl (j 2 e) = du kl + to \ A uj hl is the field strength of the Levi-Civita spin 
connection ui kl . The theory is invariant under the whole group Aut(P) meaning 
that each vector field l|47|l is an infinitesimal generator of symmetries. Field 
equations ensuing from l|49|) are, 

-LeiR k y^e ljkl = (50) 

As it is well known the theory described by (|50[1 is completely equivalent to 
Einstein's theory (provided that det(e* ) ^ 0), the difference between the two 
sets of field equations being only a matter of notation. Therefore, we should 
expect to reproduce the same expressions for the conserved quantities we derived 
in the previous section. 

The variational Lagrangian associated with the Euler-Lagrange equations 
(| 5 C)|> and the infinitesimal generator of symmetries (|47(l turns out to be: 

L'(S) = -±-e^e ijkl {elR^sel} ds (51) 

where the Lie derivative of the tetrad field has been defined by 1)48 [I. It is an easy 
calculation to show that (|51l) coincides with i|33fl so that exactly the same results 
of the previous section are recovered. Notice that, even though this is exactly 
the result we physically expected, it it nevertheless a result not so obvious from 
a mathematical viewpoint. Indeed, even if the Lagrangians (|31() and (|49|l are 
equivalent as far as their field equations contents are concerned, they describe 
instead different theories from a geometrical viewpoint, since the former theory 
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is a natural theory while the latter one is a true gauge natural theory. Therefore 
they admit two different groups of symmetries given, respectively, by Diff (M) 
and Aut(P). For this reason conserved quantities for tetrad gravity calculated 
via Noether theorem are intrinsically indeterminate from a mathematical view- 
point (see [S^I) since the components £ y of the vector field (|47|) enter into the 
Noether superpotential and they can be, a priori, fixed arbitrarily. This inde- 
terminacy does not occur in our framework since, in ^51(1. the skew-symmetric 
components £ u are contracted with the symmetric components i? Alt/ ef of the 
Ricci tensor and therefore they give rise to an identically vanishing term. 

Therefore, we expect that in any generalization of Einstein's theory involving 
also the non-symmetric part of the Ricci tensor (see, e.g. the indeter- 

minacy will still remain so that we have to somehow bypass it. This is the issue 
we shall consider in the next section. 

6 Einstein Cartan theory and the Kosmann lift 

The Einstein-Cartan (Sciama-Kibble) theory, ECKS from now on, describes a 
generalization of Einstein theory of gravitation in which the spin angular mo- 
mentum of matter plays a dynamical role. In the ECKS model the dynamical 
connection is still metric compatible but it is no longer symmetric, the torsion 
part being coupled to the spin of matter. Therefore the theory leads to devi- 
ations from General Relativity (even though the mathematical differences can 
be physically appreciated only in extreme solutions); see 

The ECKS theory can be framed in the domain of gauge natural theories; 
see ■ The configuration bundle Y for the gravitational ECKS Lagrangian is 
assumed to be the product bundle: 

F = £xC(P) — >M (52) 

where C(P) = J 1 P/50(1,3) denotes the bundle of principal connections over 
the structure bundle P (see |SH]) while £ and P have been defined in the pre- 
vious section. The bundle Y is clearly a gauge natural bundle associated with 
the structure bundle (P, M,p; 5*0(1, 3)). Fibered coordinates on Y are denoted 
by (x p ,e l ^,Al k ). The Lagrangian of the theory is the sum L = Lec + Lm 
of the gravitational ECKS Lagrangian Lec and the matter Lagrangian Lm- 
The former Lagrangian turns out to be the fibered morphism Lec '■ E x 
J 1 C(P) — > A 4 (M) locally described through the La grangian density: 

Lec = \e\ 4 F^tfA) e^e m (53) 

where F kl (j 1 A) — dA kl +A k h A A hl is the field strength of the connection (from 
now on we shall omit the constant factor l/2k in front of the Lagrangian). The 
theory is invariant under the whole group Aut(P) (provided that the matter 
Lagrangian is, in its turn, a gauge natural Lagrangian) meaning that each vector 
field H47|) is an infinitesimal generator of symmetries. Field equations ensuing 
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from 153|) through the variations of the independent fields and are, 
respectively: 

leiF^ a % kl =rt (54) 

-e^ af3 e ijkl ^ L el = rf ; (55) 

where D a e\ = 4eJ, + A 1 j Q e^, while rf and Tjy are, respectively, the energy 
momentum and the spin of matter, the explicit form of which depends on the 
matter Lagrangian. 

As it is well known, if the spin of matter is zero the theory described by l|54|l 
and Ij55(l is on-shell equivalent to Einstein's theory (provided that det(e^) 7^ 0). 
Indeed equation l(53|) admits the solution A l j a = uj l j a (j 1 e), where 0J l j a (j 1 e) 
denotes the Levi-Civita connection built out of the tetrad together with its 
derivatives. This solution inserted back into the Lagrangian l|53|) gives rise to 
the Hilbert Lagrangian (|31J) while equations l|54|) become Einstein's equations 
l|32l) . Nevertheless, even when — the equivalence holds true only on- 
shell since there does not exist any off-shell relation between the tetrad field 
and the connection. On the contrary, when the spin of matter does not vanish 
equation (|55J) becomes an algebraic equation relating the "non-metricity" of the 
connection (namely, the torsion) with the spin matter content. 

We shall now exhamine the energetic informations ensuing from the pure 
gravitational part of the theory (left hand side of equations ()54|l and (J^HJl). 
Since all calculations will be performed off-shell, the result will be not affected 
by the specific distribution of spinning matter one is taking into consideration. 
The variational Lagrangian associated with the left hand side of Euler-Lagrange 
equations (|54(1 - (|55|1 and the infinitesimal generator of symmetries Ij47(l turns out 
to be: 

L'(S) = e^e ljkl LleiF* l £ s ei + 4 D a ei£ s A k l j ds (56) 
where the Lie derivative of the tetrad field has been defined in (|48|l . while 

£ B A% = ^dpA^ + dp^A^+D^e 1 

= eF^+D?^ (^) = e fe '+<n (57) 

Notice that the Lagrangian (|56|l is a gauge natural Lagrangian by its own; see 
section |5] 

First step. The variation of l|56|) with respect to a vertical vector field X = 

Se^d/de 1 ^ + SA^^d/dA 1 ^ splits into a term < e(L')\X > which is identically 

vanishing plus a divergence term Div J-{L' : X) (see (|13f) 'l. The term T(JJ ', X) is 

A 

linear in the coefficients and their derivatives d^jD^ C 3 ■ Integrating 

by parts according to formula l|18l) and (|19|l we get (second step): 

r^(L',X) = S x ^(L,E)+d W^(L',X) (58) 
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where: 

£?(L,2) = e^e ijkl e ^ - £ Q e^)} (59) 

W?(L', X) = e^ v e ijkl ei {SA k \ L e^ + Se^ } (60) 

Therefore, the variation of conserved quantities is defined (see (|24p) as the on- 
shell integral: 

e,A) ~ \ ( WP(L>,X)ds y0 

z JdD 

- IS^^^^^e+K^ds^ (6i) 

where dD is the two dimensional boundary of a three dimensional region D. 
Notice that the components ^ of the vector field enter into the definition of 
conserved quantities. Since, up to now, no preferred condition can be math- 
ematically imposed on those components, formula (|61J) features an intrinsic 
indeterminacy. Nevertheless it was suggested in jSS] that the indeterminacy 
can be eliminated if the vertical components £^ are defined through the Kos- 
mann lift (we shall enter into details later on). The choice of the Kosmann 
lift was justified, a posteriori, in |59| from a physical viewpoint, since, in an 
apparently surprising way, it is the only lift among the possible ones which al- 
lows one to exactly reproduce the expected values for conserved quantities in 
explicit applications; see |51 I29|. Our next issue will be then to give a mathe- 
matical justification for this choice. Indeed we shall show that the Kosmann lift 
arises spontaneously when trying to restore the naturality of the gauge natural 
variational Lagrangian Ij56(l . 

To this end let us face up to the problem of calculating, through the vari- 
ational Lagrangian l|56|) . the variation of the conserved quantity associated to 
a spacetime vector field £. Notice, however, that vertical vector fields Sp = 
^(x)pij in the structure bundle (see i@3) ) are well defined objects and they de- 
scribe infinitesimal gauge transformations (and, accordingly, expression (|61|l de- 
scribes the variation of gauge charges when calculated for vertical vector fields). 
On the contrary, horizontal vector fields Sp = £ M (cc)<9 M are not globally well 
defined since they do not transform tensorially. Roughly speaking this means 
that the Lie derivatives of the dynamical fields are defined only with respect to 
vector fields in P and not with respect to spacetime vector fields. Hence, we 
have, first of all, to somehow lift the vector field £ to a vector field 3 = £ on the 
principal bundle P. If we are able to do this, the Lie derivatives £iA and £ge 

in expression (|56l) are well defined and expression (|61|l . with S = £, can be used 
to define global charges (such as energy, momentum and angular momentum). 
What we shall look for is, if it exists, a geometrically well defined lift which 
restores, off-shell, the naturality of the variational Lagrangian (|56|l . 

Notice that, as far as expression (|56|l is concerned there does not appear any 
preferred dynamical linear connection on spacetime (as it stands, expression 
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I|56|) already transforms covariantly with respect to the automorphisms of the 
principal bundle). Hence we are free to introduce, by hands and according to 

our convenience, a linear connection T built out of the dynamical fields. In this 

A 

way we can extend the SO(l, 3) covariant derivative D, which acts on internal 

r 

(Latin) indices, to a 5*0(1,3) x GL{A) covariant derivative V which acts on 
internal as well as spacctime (Greek) indices. We define the linear connection 
r 

r by requiring that V annihilates the tetrad field, i.e.: 

V 7 4 :=D, e* + rV, e* = d^ + e p + A% e£ = (62) 

These are 64 independent linear equations in the 64 variables T P AI7 . Thereby 
the coefficients of the linear connection, being uniquely determined in terms of 
the principal connection and of the tetrad field together with its derivatives, 
describe in fact a dynamical field. Notice that from Q62JI it immediately follows 
that g is parallel along T, i.e.: 

V 7 9„u = (63) 
Hence the connection turns out to be: 

rVr=7Vy(i 1 s)+^Vr ( 64 ) 

where r ) p ^{^g) denotes the Levi-Civita connection of the metric <? M „ = ^r/e^ej 
while K p M7 is the so-called contorsion tensor: 

K p /x7 — T* 3 p 7 T^~y Ty^^ 7 T p p 7 — r^[/^7] (65) 

(indices are clearly raised and lowered with the metric g^ and its inverse). In 
this step we can replace the original 24 variables A l j 7 with the 24 variables 
K p m (notice that K pfll = - A" Aip7 ). 

Equation l|62|) allows to produce a relation between the field strength F lJ ^ 
and the Riemann tensor 

R a 0p,v = rf/if 1 ^ — d v T a + r Q CTA1 r' T p v — r a au T' T i3 f j i (66) 

The relation is indeed established by the property: 

r r . . r 

Vf]cp — F j^iu^p FP pfii/^j ~t~ Y/ 7 Cp — 

— F jpu^p FO pfn/€^ (67) 
which can be easily solved in terms of the field strength as follows: 

= Wpp„ e 7 e] (68) 
Inserting (|62|l and (|68() into the variational Lagrangian l|56|) we obtain: 

L'(3) = U^(R p p - ±8* R) e P £ E el + e^e^e* e* p T>> va £ 3 A%\ ds (69) 
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where = R a mv denotes the Ricci tensor built out of the connection T. 
From (H2J it follows that: 

RputfT) = R liV {j 1 l)+ V/3 K\ v - v„ K\^K\ ? K% v -K a av K\ a (70) 

where R^vij 1 !) is the Ricci tensor of the Levi-Civita connection and V is the 
(metric) covariant derivative induced by 7. Notice that R pv is not symmetric. 
By splitting the first term in the right hand side of l|69|l into its symmetric and 
skew-symmetric parts we obtain: 

+2^/gR^e pi £ 3 ei ds (71) 

with G pp — R pp ~ l/2g pp R. The above expression suggests how to select 
a preferred lift S of the spacetime vector field £, namely, how to define the 
components £j in l|48|l in terms of the the components £ M . We simply require 
that the latter term in l|71[) does vanish, i.e: 

e^-e^ = (72) 

This algebric equation can be easily and uniquely solved in terms of the com- 
ponents £*■ yielding: 

£y = = e g e fi x d x Z a - ^egd^" (73) 

The above expression is known as the generalized Kosmann lift. We remark that 
the Kosmann lift was defined for the first time in [23 in order to establish a 
relationship between the ad hoc definition of Lie derivative of spinor fields given 
in [^Hj and the general theory of Lie derivatives on fiber bundles. We stress 
that, in our framework, the generalized Kosmann lift arises, in a completely 
independent way, when trying to restore naturality in the variational Lagrangian 
l|56l) . Namely, it is the consequence of the algebraic equation l|72ll and it is not an 
ad hoc definition. We also stress that the solution (|73|l is globally well defined. 
Indeed, the Kosmann lift K(^) of a spacetime vector field £, which in a system 
of fibered coordinates read as follows: 

K (0 = ed, + {el e^d x C - C e^e^} Pij (74) 

transforms tensorially (as can be inferred through a direct inspection) and can 
thereby globally defined. Namely, all the local expressions Ij74(l can be patched 
together to define a global vector field. 6 

6 The geometric nature of the Kosmann lift can be easily understood as follows. One starts 
with the spacetime vector field £ and considers its natutal lift to the frame bundle L(M). 
Fixing a frame e = (e*) and thereby a spacetime metric g one can construct the bundle of 
orthonormal frames SO(M, g) which is a subbundle of L(M). The Kosmann lift is nothing but 
the projection of from the tangent space of L(M) onto the tangent space of SO(M, g); 
see |41|. 
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We now analyze the consequences that the expression l|73[) has on our theory. 
Inserting (|73|l into l|48[) we have: 

£ K(C e l = \ e ' 1 ' £ H9^ (75) 

Instead if we insert l)73|l into (|57|l and we take the relations ilt)2[) and i|()8[l into 
account we get: 



£ i T\ (P (76) 



where V denotes the covariant derivatives with respect to the trasposed con 

\ _ pA 
fH — 1 n~f 



nection f A 7Al = T A M7 . Finally if we insert fT2)l. (|75j) and (J7BJ into JJTJ we end 



up with: 

(0) = y/g{G^£^ p + 2 f rT\ a £ i T\ p } ds (77) 

where: 

r\ a = T\ a - 6?TP Xp + 5 p x T? ap (78) 

is the modified torsion tensor; see |45| . Notice that the variational Lagrangian 
of our initial theory turns out to be completely natural. Notice also that (|77|l 
is nothing but the variational Lagrangian ensuing form the Einstein-Cartan 
Lagrangian density in the metric affine formalism: 

^ec = Vgg^R^V 1 ?) (79) 

where F is the metric compatible connection (I64[) . Indeed field equations ensuing 
from H79|) through a variational principle a la Palatini are nothing but the 
coefficients in front of the Lie derivatives in (|77(l (except for a sign), i.e. 



EC 



5Lec 



(80) 



7/3 



= -2^T\ a (81) 



Remark 6.1 We point out that only the symmetric components of the Ricci 
tensor enter into the Lagrangian (|79|l and in its associated variational La- 
grangian (|77|) . This fact justifies a posteriori the condition l|72|) . ■ 



Notice, however, that in order to extract the information about the varia- 
tion of conserved quantities from the variational Lagrangian 177(1 we have to 
take into account that the linear connection T is a dynamical field r^' 1 ^, .K") 
depending on the metric, its derivatives and the contorsion tensor; see (Ifi4|l . 
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In other words (when implementing the first step) the variational derivatives 
of the Lagrangian (|77)l are to be performed with respect to the independent 
fields g^v and K a f}^ (or, equivalently, with respect to g pv and T a p^). After a 
straightforward calculation we get: 



vWASw + 2g^T' 3 Xa £ i r\ /3 = Li + L 2 + L 3 (82) 



where: 



L x = s /gG^(fg)£ £ g^ (83) 
L2 = £t;{-V99^(K p * P K%„~K a au K^ a )} (84) 
L 3 = d p {-2^T^£ i9fw )} (85) 

The potential for L\ has already been computed in section^and is given (apart 
from a factor 2k) by formula l|42|l. Moreover there is no contribution to the 
potential from L 2 since L 2 affects only the term 58 in i|23|) . Indeed we have 
(first step): 

5L 2 = d p {-^5 {^ v {K a aa K\ v - K a av K^ a ))} (86) 
hence, see (fTT?|) : 

F p (L 2 ,6g, 5K) = -£'<S (y/gg^ \K a aa K" ^ - K a au K^ a )) (87) 

The latter expression is linear in £ so that, in the second step, no further inte- 
gration with respect to £ can be made, namely no divergence term entering into 
the total potential comes out from it. 

The potential receives instead a contribution from the term L3. To obtain 
it we have to consider the part enclosed between the braces in Q85J1. develop 
explicitly the Lie derivatives and then apply fTTHl with h + r—1 = 1. The result 
turns out to be 

U 0a (L 3l 6g,6T) = 6 [iy/g^^ + W} (88) 

Hence the total potential for the variational Lagrangian (|77|l is given by the sum 
of the purely metric potential (|42|l plus 188(1 . namely: 

+ s\2^g(T^+T p al3 )C p } (89) 

where u£„ = 7^„ — 5? 7^ CT J We remark that, if the spin contribution of the 
matter is zero, it holds true that T p \ a ~ 0. Hence, expression ()89fl reproduces 
on-shell exactly expression l|42ll . When the spin contribution of matter is not 
vanishing it gives to conserved quantities a geometric contribution (i.e. relative 

7 The same result 1891 could be also achieved, in a quite more involved way, by inserting 
the expression 1731 of the Kosmann lift directly into H6UI 
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to the pure gravitational part of the theory) which is produced by the torsion 
terms in (|89|) . Notice however that the equation l|55|l is an algebric relation be- 
tween torsion and spin matter meaning that, outside the distribution of matter, 
there is no torsion, i.e. torsion does not propagate (see @5J). The only way the 
spin of matter affects conserved quantities outside the distribution of matter 
(where collapses into is through its influence on the metric tensor. 

7 Conclusions and perspectives 

We have described a theory which algorithmically defines the variation of con- 
served quantities in the realm of gauge natural theories. The formalism is de- 
veloped through a two steps procedure starting from the variational Lagrangian 
and by making use of techniques related to the Calculus of Variations in the 
jet bundle framework. We have tested the viability of the results the theory 
predicts by applying them to gravitational theories. In the classical Einstein's 
metric formulation we reproduced results already found elsewhere. Moreover, 
when dealing with ECKS theory, viewed as a gauge natural generalization of 
gravity, a preferred lift of spacetime vector fields, namely the Kosmann lift, is 
automatically selected by the formalism itself in order to restore naturality in 
the gauge natural framework. We clearly welcome this additional unlooked-for 
result which has turned up from the theory since it seems to endow the Kos- 
mann lift, that in literature was previously justified only on physical grounds 
(see PJIJniESl), also with a good mathematical reason of existence. 

Moreover we point out that the left hand side of equations (|51J) and 1 (55) 1 
can be easily generalized to encompass any spacetime dimension D, D > 2. In 
particular, their generalization to odd dimensions reproduces the field equations 
content of Chern-Simon theories of the group SO(2, D — 1) or SO(l, D), see |T8) . 
The case D — 3 has been analyzed in [3] following the spirit of the formalism 
presented here. We believe that it will be worth investigating also the case 
D = 5 (as well as any larger dimension), mainly because of its relationship with 
the Gauss-Bonnet generalization of Einstein's gravity, a generalization which 
seems to be relevant in effective theories ensuing from the low-energy limit of 
string theory (see 03] ). 

We also stress that our recipe is well suited to be generalized to a more gen- 
eral class of symmetries other than those induced just by projectable vector fieds 
on the structure bundle. Indeed the fundamental requirement we made on the 
Lie derivatives of the fields is that they can be expanded as a linear combination 
of the independent parameters £^ and £ A together with their derivatives up to 
a fixed (finite) order. This is in fact the assumption which makes the formal- 
ism algorithmically well-defined (see expression © and (|I7Jl ). No restriction 
has been instead imposed on the coefficients of the aforementioned expansion 
which, in principle, can depend not only on the field together with their first 
order derivatives, but also on an arbitrary (but finite) number of derivatives. 
Namely, symmetries with respect to generalized vector fields (see j2H E2] are 
not out of the scope of the formalism here developed. Thereby we guess that the 



2G 



formalism is in good position to deal with the so-called generalized symmetries, 
e.g. supersymmetries and BRST transformations; see |32l I46| . 
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